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NON-EQUILIBRIUM ISOTHERMAL TRANSFORMATIONS IN A 
TEMPERATURE GRADIENT FROM A MICROSCOPIC DYNAMICS 

VIVIANA LETIZIA AND STEFANO OLLA 


Abstract. We consider a chain of anharmonic oscillators immersed in a heat bath 
with a temperature gradient and a time-varying tension applied to one end of the 
chain while the other side is fixed to a point. We prove that under diffusive space-time 
rescaling the volume strain distribution of the chain evolves following a non-linear 
diffusive equation. The stationary states of the dynamics are of non-equilibrium and 
have a positive entropy production, so the classical relative entropy methods cannot 
be used. We develop new estimates based on entropic hypocoercivity, that allow to 
control the distribution of the position configurations of the chain. The macroscopic 
limit can be used to model isothermal thermodynamic transformations between non¬ 
equilibrium stationary states. 


1. Introduction 

Macroscopic isothermal thermodynamic transformations can be modeled microscop¬ 
ically by putting a system in contact with Langevin heat bath at a given temperature 
/3-1. In m a chain of n anharmonic oscillators is immersed in a heat bath of Langevin 
thermostats acting independently on each particle. Macroscopically equivalent isother¬ 
mal dynamics is obtained by elastic collisions with an external gas of independent 
particles with maxwellian random velocities with variance /3~^. The effect is to quickly 
renew the velocities distribution of the particles, so that at any given time it is very 
close to a maxwellian at given temperature. The chain is pinned only on one side, while 
at the opposite site a force (tension) r is acting. The equilibrium distribution is char¬ 
acterized by the two control parameters /3“^,r (temperature and tension). The total 
length and the energy of the system in equilibrium are in general non-linear functions 
of these parameters given by the standard thermodynamic relations. 

By changing the tension r applied to the system, a new equilibrium state, with the 
same temperature will be eventually reached. For large n, while the heat bath equi¬ 
librates the velocities at the corresponding temperature at time of order 1, the system 
converges to this global equilibrium length at a time scale of order n^t. In it is proven 
that the length stretch of the system evolves in a diffusive space-time scale, i.e. after 
a scaling limit the empirical distribution of the interparticle distances converges to the 
solution of a non-linear diffusive equation governed by the local tension. Consequently 
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this diffusive equation describes the non-reversible isothermal thermodynamic transfor¬ 
mation from one equilibrium to another with a different tension. By a further rescaling 
of the time dependence of the changing tension, a so called quasi-static or reversible 
isothermal transformation is obtained. Corresponding Clausius equalities/inequalities 
relating work done and change in free energy can be proven. 

The results of [9] summarized above concern isothermal transformations from an 
equilibrium state to another, by changing the applied tension. In this article we are 
interested in transformations between non-equilibrium stationary states. We now con¬ 
sider the chain of oscillators immersed in a heat bath with a maeroseopic gradient of 
temperature: each particle is in contact with thermostats at a different temperature. 
These temperatures slowly change from a particle to the neighboring one. A tension 
r is again applied to the chain. In the stationary state, that is now characterized by 
the tension r and the profile of temperatures there is a continuous flow 

of energy through the chain from the hot thermostats to the cold ones. Unlike the 
equilibrium case, the probability distribution of the configurations of the chain in the 
stationary state cannot be computed explicitly. 

By changing the applied tension we can obtain transitions from a non-equilibrium 
stationary state to another, that will happen in a diffusive space-time scale as in the 
equilibrium case. The main result in the present article is that these transformations 
are again governed by a diffusive equation that takes into account the local temperature 
profile. The free energy can be computed according to the local equilibrium rule and 
its changes during the transformation satisfy the Clausius inequality with respect to 
the work done. This provides a mathematically precise example for understanding 
non-equilibrium thermodynamics from microscopic dynamics. 

The results in [9] where obtained by using the relative entropy method, first devel¬ 
oped by H.T.Yau in m for the Ginzburg-Landau dynamics, which is just the over¬ 
damped version of the bulk dynamics of the oscillators chain. The relative entropy 
method is very powerful and flexible, and was already applied to interacting Ornstein- 
Uhlenbeck particles in the PhD thesis of Tremoulet m as well as many other cases, 
in particular in the hyperbolic scaling limit for Euler equation in the smooth regime 
[lam- This method consists in looking at the time evolution of the relative entropy of 
the distribution of the particle with respect to the local Gibbs measure parametrized 
by the non-constant tension profile corresponding to the solution of the macroscopic 
diffusion equation. The point of the method is in proving that the time derivative of 
such relative entropy is small, so that the relative entropy itself remains small with 
respect to the size of the system and local equilibrium, in a weak but sufficient form, 
propagates in time. In the particular applications to interacting Ornstein-Uhlenbeck 
particles [Hi, the local Gibbs measure needs to be corrected by a small recentering 
of the damped velocities due to the local gradient of the tension. 

The relative entropy method seems to fail when the stationary measures are not the 
equilibrium Gibbs measure, like in the present case. The reason is that when taking the 
time derivative of the relative entropy mentioned above, a large term, proportional to 
the gradient of the temperature, appears. This term is related to the entropy produe- 
tion of the stationary measure. Consequently we could not apply the relative entropy 
method to the present problem. 

A previous method was developed by Guo, Papanicolaou and Varadhan in [6] for over¬ 
damped dynamics. In this approach the main step in closing the macroscopic equation 
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is the direct comparison of the coarse grained empirical density in the microscopic and 
macroscopic space scale. They obtain first a bound of the Dirichlet form (more precisely 
called Fisher information) from the time derivative of the relative entropy with respect 
to the equilibrium stationary measures. This bound implies that the system is close to 
equilibrium on a local microscopic scale, and that the density on a large microscopic 
interval is close to the density in a small macroscopic interval (the so called one and 
two block estimates, see [7] chapter 5). 

In the over-damped dynamics considered in [6], the Dirichlet form appearing in 
the time derivative of the relative entropy controls the gradients of the probability 
distributions with respects to the position of the particles. In the damped models, 
the Dirichlet form appearing in the time derivative of the relative entropy controls 
only the gradients on the velocities of the probability distribution of the particles. In 
order to deal with damped models a different approach for comparing densities on the 
different scales was developed in [12], after the over-damped case in m, based on 
Young measures. Unfortunately this approach requires a control of higher moments of 
the density that are difficult to prove for lattice models. Consequently we could not 
apply this method either in the present situation. 

The main mathematical novelty in the present article is the use of entropic hypoco- 
ercivity, inspired by [l6|. We introduce a Fisher information form associated to the 
vector fields {^p^ -|- defined by (|2.27p . By computing the time derivative 

of this Fisher information form on the distribution at time t of the configurations, we 
obtain a uniform bound In < Cn~^. This implies that, at the macroscopic diffusive 
time scale, velocity gradients of the distribution are very close to positions gradients. 
This allows to obtain a bound on the Fisher information on the positions from the 
bound on the Fisher information on the velocities. At this point we are essentially with 
the same information as in the over-damped model, and we proceed as in |6]. Observe 
that the Fisher information we introduce in (I2.27h is more specific and a bit different 
than the distorted Fisher information used by Villani in [T6|, in particular is more 
degenerate. On the other hand the calculations, that are contained in appendix D are 
less miraculous than in m, and they are stable enough to control the effect of the 
boundary tension and of the gradient of temperature. This also suggests that entropic 
hypocoercivity seems to be the right tool in order to obtain explicit estimates uniform 
in the dimension of the system. 

Adiabatic thermodynamic transformations are certainly more difficult to be obtained 
from microscopic dynamics, for some preliminary results see [laSKIKlo]. Equilibrium 
fluctuations for the dynamics with constant temperature can be treated as in HH. The 
fluctuations in the case with a gradient of temperature are non-equilibrium fluctuation, 
and we believe that can be treated with the techniques of the present article together 
with those developed in the over-damped case in [5|. 

Large deviations for the stationary measure also require some further mathematical 
investigations, but we conjecture that the corresponding quasi-potential functional (|2|) 
is given by the free energy associated to the local Gibbs measure, without any non-local 
terms, unlike the case of the simple exclusion process. 

The article is structured in the following way. In isection 2l we dehne the dynamics and 
we state the main result (Theorem I 2.11) . In [section 3) we discuss the consequences for 
the thermodynamic transformations from a stationary state to another, the Clausius 
inequality and the quasi-static limit. In [section 4l are obtained the bounds on the 
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entropy and the various Fisher informations needed in the proof of the hydrodynamic 
limit. In [section ^ we show that any limit point of the distribution of the empirical 
density on strain of the volume is concentrated in the weak solutions of the macroscopic 
diffusion equation. The compactness, regularity and uniqueness of the corresponding 
weak solution, necessary to conclude the proof, are proven in the first three appendices. 
Appendix D contains the calculations and estimates for the time derivative of the Fisher 
information 7„. 


2. The dynamics and the results 

We consider a chain of n coupled oscillators in one dimension. Each particle has 
the same mass, equal to one. The conhguration in the phase space is described by 
i = 1,..., n} G M^"'. The interaction between two particles i and i—1 is described 
by the potential energy V{qi — qi-i) of an anharmonic spring. The chain is attached on 
the left to a fixed point, so we set go = 0,po = 0. We call {rj = qi — qi-i,i = 1,..., n} 
the interparticle distance. 


We assume V to be a positive smooth function, that satisfy the following assumptions: 

i) 


v V{r) 
lim —r—^ = oc, 

(2.1) 

|r|— )-oo \t\ 


ii) there exists a constant C 2 > 0 such that: 


sup|E"(r)| <C2, 

r 

(2.2) 

iii) there exists a constant (7i > 0 such that: 


V'{rf <Ci[l + V{r))). 

(2.3) 


In particular these conditions imply |F'(r)| < Cq + (721^1 for some constant Cq. Notice 
that this conditions allows potentials that may grow like V{r) ~ |r|“ for large r, with 
1 < a < 2. 

Energy is defined by the following Hamiltonian function: 

^ / 2 \ 

Y + (2.4) 

i=l ^ ^ 

The particle dynamics is subject to an interaction with an environment given by 
Langevin heat bath at different temperatures We choose /3j as slowly varying on 
a macroscopic scale, i.e. /3i = /3(z/n) for a given smooth strictly positive function j3{x), 
X G [0,1] such that infj^gjoq] /3(i/) > /3_ > 0. 

The equations of motion are given by 

'dri{t) = 'n?{pi{t) -pi_i{t))dt 

< dpi{t) = 'n?{V'{ri+i{t)) - V'{ri{t))) dt - n'^'ypi{t)dt + n^J^dwi{t), i = 1, ..,N - 1 
dpn{t) = - V'{rn{t))) dt - n^^pn{t) dt + n^J^dwn{t). 

(2.5) 

Here {wi{t)}i are n-independent Wiener processes, 7 > 0 is the coupling parameter 
with the Langevin thermostats. The time is rescaled according to the diffusive space- 
time scaling, i.e. t is the macroscopic time. The tension f = f{t) changes at the 
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macroscopic time scale (i.e. very slowly in the microscopic time scale). The generator 
of the diffusion is given by 

/:;(*):= (2.6) 

where is the Liouville generator 

n n—1 

K = + Y.^V\n+^) - V'{n))dp^ + (f- y'(rO)ap„ (2.7) 

i=l i=l 

while Sn is the operator 

n 

Sn = Y^ - pidp^) ( 2 . 8 ) 

i=l 


2.1. Gibbs measures. For f{t) = r constant, and f3i = [3 homogeneous, the sys¬ 
tem has a unique invariant probability measure given by a product of invariant Gibbs 
measures //” ^: 

n 

dPr,p = n dndpi (2.9) 

i=l 

where £i is the energy of the particle i: 

2 

^* = Y + ^(^*)- (2-10) 

The function Q{T,f3) is the Gibbs potential defined as: 


Q{t,15) = log 


^ 271 / 3-1 J 


( 2 . 11 ) 


Notice that, thanks to condition (j2.1h . Q{t,I3) is finite for any r G M and any /3 > 0. 
Furthermore it is strictly convex in r. 

The free energy of the equilibrium state (r, (3) is given by the Legendre transform of 

r'e(r,/3): 

J-(r, /3) = sup{rr - /3-^g(r, /?)} (2.12) 

r 

The corresponding convex conjugate variables are the equilibrium average length 

t(G Id) = l3~^drG{T, (3) (2.13) 

and the tension 

T(r,/3) = a^J’(r,/3). (2.14) 


Observe that 


%,,h]=t(r,/3), E^.^[G'(r,)]=r. 


(2.15) 
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2.2. The hydrodynamic limit. We assume that for a given initial profile ro(x) the 
initial probability distribution satisfies: 

— } G(i/n)ri(0) —> / G{x)ro(x)dx in probability (2-16) 

n n^co /„ 

i=l 

for any continuous test function G E Co([0,1]). We expect that this same convergence 
happens at the macroscopic time t: 


1 ” G 

— G{i/n)ri{t )—/ G(x)r(x,t)dx 

^ i=i 


(2.17) 


where r{x,t) satisfies the following diffusive equation 

dtr{x,t) = ^d‘^T{r{x,t),l3{x)) for x E [0,1] 

dxT{r{t,x),(3{x))\x=o = 0, T{r{t,x),(3{x))\x=i = f{t), t>0 (2.18) 

r(0, x) = ro(x), x E [0,1] 

We say that r{x,t) is a weak solution of (12.181) if for any smooth function G{x) on 
[0,1] such that G(l) = 0 and G'(0) = 0 we have 

f G{x) {r{x,t) — ro{x)) dx = f ds f G"{x)T{r{x, s), (3{x))dx — G'{l)f{s) 

Jo Jo Uo 

(2.19) 

In appendix C we prove that the weak solution is unique in the class of functions such 
that: 

I ds {dxT{r{x, s), (3{x)))'^ dx <+ 00 . (2.20) 

Jo Jo 

Let the inhomogeneous Gibbs measure 


■i"?.=n 


n —Rf- 

, _^ g 


2 = 1 


% 


( 2 . 21 ) 


Observe that this is not the stationary measure for the dynamics defined by (USD 
and (j2.6jl for r = 0. 

Let the density, with respect to , of the probability distribution of the system 
at time t, i.e. the solution of 

dtf? = 

where is the adjoint of with respect to v'^ , i.e. explicitly 


( 2 . 22 ) 


n—1 


^ -n'Y^ VnlJ{i/n)piV'{ri+i) + n^/3(l)p„r + n^75„, (2.23) 


2 = 1 


where 


V„/3(z/n) = n (/3 


i + 1 


n 


-/3 


n 


i = 1,..., n — 1. 


Define the relative entropy of ftdv^ with respect to dv'p as: 

Hn{t) = I log f^du^. 


(2.24) 


(2.25) 
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We assume that the initial density satisfy the bound 

Hn{0) < Cn. (2.26) 

We also need some regularity of Jq: define the hypercoercive Fisher information func¬ 
tional: 

E A' f (2.27) 

i=i 

where = dr^ — Sri+u ^ = 1,..., n — 1, and := We assume that 

4(0) < Kn (2.28) 


with Kji growing less than exponentially in n. We will show in Appendix D that for 
any t > 0 we have /n(t) < Cn~^. 

Furthermore we assume that 


lim 

n—)-oo 


I 1 GW'-oW* 


fSdi,f. = 0 


(2.29) 


for any continuous test function G G Co([0, 1]). 


Theorem 2.1. Assume that the starting initial distribution satisfy the above conditions. 
Then 


lim 

n^oo 



f?diyg. = 0, 


(2.30) 


where r{x,t) is the unique weak solution of \2.18\) satisfying (|2.2Up . 


Furthermore a local equilibrium result is valid in the following sense: consider a local 
function (^(r,p) such that for some positive finite constants Ci,C 2 we have the bound 


\(l){r,p)\ +V{ri))°‘+ C 2 , a<l (2.31) 

where A^p is the local support of 4>. Let kp, the length of A<^, and let 6i4> be the shifted 
function, well dehned for < z < n — k^p, and define 


= ^U.ir,p),p (<^) 


(2.32) 


Corollary 2.2. 


lim 

n—>-oo 




n 


E G 

i=k^+l 


9i(f){r,p) - [ G{x)^{r{x,t),fi{x))dx 
nj Jo 


f'tdug. = 0, (2.33) 


3. Non-equilibrium thermodynamics 

We collect in this section some interesting consequences of the main theorem for 
the non-equilibrium thermodynamics of this system. All statements contained in this 
section can be proven rigorously, except for one that will require more investigation 
in the future. The aim is to build a non equilibrium thermodynamics in the spirit of 
[HI [2] . The equilibrium version of these results has been already proven in |9] . 
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As we already mentioned, stationary states of our dynamics are not given by Gibbs 
measures if a gradient in the temperature profile is present, but they are still char¬ 
acterized by the tension f applied. We denote these stationary distributions as non¬ 
equilibrium stationary states (NESS). Let us denote the density of the stationary 
distribution with respect to . 

It is easy to see that 

j y'{ri)fgg.^vp, =T, i = l,...,n. (3.1) 

In fact, since f Pif^g = 0 and 

n~‘^Clpi = V'{ri+i) - V'{ri) - 'ypi, i = 1,..., n - 1, 

n~‘^Clpn = T - V'{rn) - IPn, 


we have 

0 = j{V'{r,+^)-V'{n))fl^,vp,= j{T-V'{rn))fZ,r^p.. 

By the main theorem l2.ll there exists a stationary profile of stretch rgg^riu) = I3{y)) 

(defined by (I2.13p l such that for any continuous test function G: 


lim 

n—>-CxD 


/ 


i ~) ~ / G{x)rgg^r{x)dx 


n 


i=l 


n 


10 


fss,Tdi^ft. = 0 > 


(3.2) 


In order to study the transition from one stationary state to another with different 
tension, we start the system at time 0 with a stationary state with tension tq, and 
we change tension with time, setting f(f) = ri for t > ti. The distribution of the 
system will eventually converge to a stationary state with tension ti. Let r{x,t) be 
the solution of the macroscopic equation (|2.19p starting with ro(x) = rss^T-o(a;)- Clearly 
r{x, t) —5- ri(x) = rss,Ti (x), as t ^ oo. 


3.1. Excess Heat. The (normalized) total internal energy of the system is defined by 

1^/2 \ 

(I+ ''(’■•)) P-3) 

It evolves as: 

Unit) -UniO)=Wn{t)+Quit) 

where 

'^nit) = [ fis)npnis)ds = [ r(s)^^^^^ 

Jo Jo n 

is the (normalized) work done by the force r(s) up to time t, while 

Qn{t)='yn'^ ds{p){s) - I Pjis)dwiis). (3.4) 

i=i i=i 


is the total flux of energy between the system and the heat bath (divided by n). As a 
consequence of theorem I 2.11 we have that 
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where C{t) = r(x, t)dx, the total macroscopic length at time t. While for the energy 
difference we expect that 

lim (Unit) - UniO)) = [ [u{T{r{x,t),f5{x)),(3{x))-u{To,P{x))]dx (3.5) 

n^oo Jq 

where u{t,I3) is the average energy for i-e. 

u{t,/3) = I = ^ + / 

with G{t,I3) = log/Unfortunately (j3.5jl does not follow from (|2.33p . 
since (j2.3ip is not satished. Consequently at the moment we do not have a rigorous 
proof of (j3.5jl . In the constant temperature profile case, treated in [9], this limit can be 
computed rigorously thanks to the use on the relative entropy method m that gives 
a better control on the local distribution of the energy. 

Since T{r{x,t), I3{x)) —>■ ti as t —>■ oo, it follows that 

u{T{r{x,t),P{x)),/3{x)) u{n,j3{x)) 

and the energy change will become 

/*! /*+oo 

/ I3{x)) — u{to,/ 3{x))) dx = / f{s)dC{s)ds + Q = W + Q (3.6) 

Jo Jo 

where Q is the limit of (13.4p . which is called excess heat. So equation (13.6|] is the 
expression of the first principle of thermodynamics in this isothermal transformation 
between non-equilibrium stationary states. Here isothermal means that the profile of 
temperature does not change in time during the transformation. 


3.2. Free energy. Define the free energy associated to the macroscopic profile r{x,t): 

iF{t) = [ F{r{x,t), fJ{x))dx. (3.7) 

Jo 

Correspondingly the free energy associated to the macroscopic stationary state is: 

^ss{t) = f iF{rss,Tix),IJ{x))dx (3.8) 

Jo 


/o 

A straightforward calculation using (I2.19P gives 

rt 


T{t) - Tssiro) = W{t) - j ^ f ds [ {dxT{r{x,s),^{x))f dx 

Jo Jo 


and after the time limit t —)• oo 

^ss 


^ /*+oo 1*1 

(n) --^ss('ro) = W - 7"W dt {dxT{r{x,t),f3{x)))‘^ dx 

Jo Jo 


(3.9) 


(3.10) 


< W 


i.e. Clausius inequality for NESS. Notice that in the case /Jj constant, this is just the 
usual Clausius inequality (see m)- 
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3.3. Quasi-static limit and reversible transformations. The thermodynamic trans¬ 
formation obtained above from the stationary state at tension tq to the one at tension 
Ti is an irreversible transformation, where the work done on the system by the external 
force is strictly bigger than the change in free energy. 

In thermodynamics the quasi-static transformations are (vaguely) defined as those 
processes where changes are so slow such that the system is in equilibrium at each 
instant of time. In the spirit of [3] and [9], these quasi static transformations are 
precisely defined as a limiting process by rescaling the time dependence of the driving 
tension r by a small parameter e, i.e. by choosing f{et). Of course the right time scale 
at which the evolution appears is £~^t and the rescaled solution f’^(x,t) = r{x,£~^t) 
satisfy the equation 

dtf^{x, t) = ^d^T(r^{x, t),P{x)) for x € [ 0 , 1 ] 

< dxT{f'^{t,x),l3{x))\x=o = 0, T{r^{t,x),/3{x))\x=i = f{t), t>0 (3.11) 

_ t(F(0,x),/3(x)) = To, xG[ 0 , 1 ] 

By repeating the argument above, equation (I3.10|) became: 

^ ^ 1 ^- 1-00 ^1 

j;s(n)--^ 55 (^ 0 ) = -/ dt {dxT{f^x,t),l3{x))f dx (3.12) 

€7 Jo Jo 

By the same argument used in [9] for (3 constant, it can be proven that the last term on 
the right hand side of (j3.12p converges to 0 as e —)• 0, and that t(P(x, t),/3(x)) —)■ f{t) 
for almost any x € [0,1] and t > 0. Consequently in the quasi-static limit we have the 
Clausius equality 

-^.s(ri) - .F,,(ro) = W 

This implies the following equality for the heat in the quasi-static limit: 

Q= / f3-\x){S{rss{x , Tl), Ugs (x, Tl) ) S(t ssix iXq) ^Ugsix dx (3.13) 

70 

analogous of the equilibrium equality Q = TAS. 

In [ 8 ] a direct quasi-static limit is obtained form the microscopic dynamics without 
passing through the macroscopic equation ()2.19p . by choosing a driving tension f that 
changes at a slower time scale. 

4. Entropy and hypercoercive bounds 

In this section we prove the bounds on the relative entropy and the different Fisher 
informations that we need in the proof of the hydrodynamic limit in section [section 5l 
These bounds provide a quantitative information on the closeness of the local distribu¬ 
tions of the particles to some equilibrium measure. 

In order to shorten formulas, we introduce here some vectorial notation. Given two 
vectors u = {ui,... ,Un),v = {vi,... ,Vn), dehne 

n n—1 

uQv = '^(3~^UiVi, uQv = '^l3~^UiVi, \u\q = uQu, |u|| = uQu. 
i=l i=l 

We also use the notations 

= (9pi,..., dpj d* = {d*^,..., dpj, dp. = /3iPi - dp, 

dg = {dg„. . .,dgj, dg, = 5^. - . 


(4.1) 
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Observe that with this notations we can write 

Sn = -dp Qdp, Al=p-dq- dqV ■ dp + rdp^ 


(4.2) 


where V = Yli ^(^0 the • denotes the usual scalar product in R”. Then we define 

the following Fisher informations forms on a probability density distribution (with 


respect to vp.)'. 


Kif) = J 


Kif) = J 

cy-. 

InU) = J 

f \dpf + dqf\- 

' f 


\dpf\ 




(4.3) 


dup.=VPM)+Kif) + ^ 


dqfQdpf 

f 


dvji. > 0 


Proposition 4.1. Let the solution of the forward equation (12.221) . Then there exist 
a constant C such that 


Hn{t) < Cn, 


/o n 


fK{C)ds<-. 

Jo n 


(4.4) 


Proof. Taking the time derivative of the entropy we obtain: 

4h„(«)= 


(4.5) 


So that, using (I2.23p . we have 

j^Hnit) = I /r/:;«log/rdu^. = I -7n2lP^(/r) 

n—1 „ „ 

= -n'^Vn/3{i/n) j V {ri+i)pifjfdvp, + (3nf{t) J Pnffdu/s,-jn^VP{fi 


2=1 


t J 

(4.6) 


Recall that Qn = YJa=i then the time integral of the second term on the RHS of (|4.6p 
gives 


n 


’(3n J ds r(s) J PnfsdJ^is. = J ds t{s) j C^'^'^qnfsdvp. 


(4.7) 


= /3nr(t) J qnf?dv 0 , - /3nr(0) J qnfodvp, ~ J 

By the entropy inequality, for any ai > 0, using the first of the conditions (|2.ip . 

[ iQnlfsdiy/s. < — log [ + —Hnis) < — log [ T\e'^^^''^^di^p. + —Hn{s) 

J ai J ai ai J oi 

2=1 

< — y log / (e“i"* + dng + —77„(s) 

ai ^ J ' ^ ai 

1=1 

1 1 1 

= — y {G{ai,Pi) + Q(-ai,(3i) - 2^(0,/3i)) H- Hn{s) < nC{ai,/3.) H- Hn{s) 


2=1 


(4.8) 
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We apply (j4.8p to the three terms of the RHS of (14.7p . So after this time integration 
we can estimate, for any oi > 0 , 


n'/3(l) 


ds T{t) / Pnftdvii. 


< 


mKf 


ai 




M +r 

40 


Hn{s)ds 


(4.9) 


+n (2 + t)/3(l)i^^C(ai,/3.) 

where = sup^>o (l'^('S)l + l^'('S)l). 

By integration by part and Schwarz inequality, for any 02 > 0 we have 


n—1 


n'^Vn^{i/n) / V'{ri+i)pif^dup, 


2=1 


n 


n—l 

E 


Vn/3(i/n) 


< 


- n—l 

-E 




2a2 ^ 


2=1 


A 


I V'{r,+,) 


V’{ri+i)dpJl^di^fS, 


"frill's. + '^m/r) 


By our assumptions on I3{-) and assumption (12.3p on V, we have that for some 
constant C/j. > 0 depending on /3(-) and V, 

n—1 /ry Q(' / ^^2 ^^ ^ 

E y'{n+if < Cp. Y, V'{n+if < CpCi Y + 1) (4-10) 

i=l i=l i=l 

By the entropy inequality, for any 5 such that 0 < 5 < infy/3(?/), there exists a finite 
constant depending on V,5 and j3{-) such that: 


Y [ V{ri)f^dup. <]log f + 

i=l 


1 ^ 1 1 

J (^(0) A - (5) - ^(0,/3i)) + -Hn{t) < + -Hn{t) 


(4.11) 

At this point we have obtained the following inequality, for some constant C not 
depending on n, 

H^{t) - Hn{0) < -r? (7 - f VPifs)ds + H^is)ds 


{Hn{t) + HniO)) + nc(ai, 02 ,5,r,/?.) 


ai 


consequently, choosing 02 = 7 and oi = 2/3{l)Kf, we have 


Hjiit) < 327, 


,(0) + C" fHn{s)d. 

Jo 


s + cn-n‘^j VP{f^)ds 


(4.12) 


(4.13) 


where C and c are constants independent of n. Given the initial bound on 27^(0) < cn, 
by Gronwall inequality we have for some c" independent on n: 

Hn{t) < d'e^^n. 

Inserting this in p4.13p we obtain, for some C independent of n, 


7 f ‘ K ( ads <^ 


n 


(4.14) 

(4.15) 
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The bound (|4.15p gives only informations about the distribution of the velocities, 
but we actually need a corresponding bound of the distribution of the positions. 

In appendix D we prove that, as a consequence of (I4.15p . we have 


Consequently 


. / X C 

In{t) < — 

n 


Vt > 0. 


KUt 


ujn - f’S(/,") - 2 / di-0- 

<-- KifT) + + ipK/,") 

n 2 


that gives 

2C 

p;(/r) < 2VPMn + — 

n 

Since we have already the bound (I4.15p . (14.4p follows. 


(4.16) 


□ 


5. Characterization of the limit points 
Define the empirical measure 

1 

Titidx) := - 

n ' 

i=\ 

and we use the notation, for a given smooth function G : [0,1] —)• M, 

W.G):=iX;o (-)’■■(*) 

n \n 

1=1 ' ^ 

Computing the time derivative we have: 

«,G)-{<,G)= (51) 

Jo ™ \^/ 

Since 

a^^hi = n^{pi-pi_i), i = l,...,n, po = 0, 
after performing a summation by parts, we obtain 

n-l , . X 

£:W(7rr,G) = -^V„G(^^Jpi(t) + np„(f)G(l). (5.2) 

where VnG is dehned by (|2.24p . We dehne also 



i = 2,... ,n. 
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Now observe that 

n—l 


pT{t) 


\y^^nG(-]pi--pnG{l) 

\nj n 


2 = 1 
n—l 


n—l 


2=1 


7 VnG ( - ) Pi + 7 np„G(l) 


2 = 1 
n—l 

= -lY.^nG 
2=1 
^ n—l 


+ ^ V„G f - j (y'(ri+i) - V'{n)) - nG (1) (f (t) - F'K)) 


n 


Pi + ^npnG{l) 


i=2 


+ -y] KVnG - V'(ri+i) + VnG 
n \n 


n — l 


n 


l"'(rn)-V„G(-)y'(ri) 


- nG{l) (r(t) - y'(r„)) 


(5.3) 

Recall that, by the weak formulation of the macroscopic equation, cf. (j2.19p . it is 
enough to consider test functions G such that G(l) = 0 and G'(0) = 0. This takes care 
of the last term on the RHS of the above expression and in (I5.2I) . and putting these 
two expression together and dividing by 7 , we obtain 

n—l 


/:;W( 7 r-,G) = i-^(-V;V„)G ( 1 ) R'(r,+i) - 7 -V„G 

yn ^ V n / 




n — l 


n 


V'{rn) 


n—l 


+ 7 -'VnG - R'(ri) + VnG - 

\n) 7 n^ ^ \ n 


2=1 


(5.4) 


Pi 


R is easy to show, by using the entropy inequality, that the last two terms are negligible. 
In fact, since G'(0) = 0 we have that |VnG ( 7 ) | < Gcn”^. Furthermore 


/ 


g«T'WI < +00 Va > 0 . 


Then, using the entropy inequality we have for any a > 0: 
y-'VnG R'(ri) 

J n^a 


/ 


nja 


f^di^p. < — / \V'{r^)\f^dup, 

nj J 

G{a) ^ G' 

n a 


(5.5) 


that goes to 0 after taking the limit as n —>■ 00 then a ^ 00 . About the last term of 
the RHS in (j5.4p . after time integration we have to estimate 

n—l 


1 

7n" 


E 

2=1 


VnG 


n 


\pi\fyi'i3. 


for s = 0,t. By similar use of the entropy inequality it follows that also this term 
disappear when n —)• 00 . 

To deal with the second term of the RHS of (j5.4p . we need the following lemma: 


Lemma 5.1. 


lim E 

n—>-oo 



(R'(rn(s)) - r(s)) ds 


= 0 


(5.6) 
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Proof. Observe that 

V'{rn) - f{s) = - - IPn = - ^C^^'^Kpn + IQu)- (5.7) 

Then after time integration: 

[ {V'{rn{s)) - f{s)) ds= ^ (Pn(0) - Pn{t)) - ^{qn{t) “ ^^(0)) + 

7 q n 

It is easy to show that, using similar estimate as (SH) and (j4.8p . the expectation 
of the absolute value of the right hand side of the above expression converges to 0 as 
n oo. □ 


It follows that 


lim E 

n—>-oo 


V„,G 


n — 1 


n 


^'(^n(s)) - G'(l)r(s) ) d,s 


= 0 . 


(5.8) 


We are finally left to deal with the first term of the RHS of (j5.4p . We will proceed 
as in [6]. For any e > 0 define 

2 ne + l T] 

\j-i\<n£ 


Vi.E = 


ne <i < n(l — s). 


(5.9) 


We first prove that the boundary terms are negligible: 

Lemma 5.2. 


lim lim 

£—>0 n—>co 



/ [ne] [n-1] 

E+ E I (-v;v„)G ( 1 ) U(r„) 


yn 


\i=2 j=[n(l-£)]+l^ 


/” dvp, ds = 0 
(5.10) 


Proof. For simplicity of notation let us estimate just one side. Since our conditions on 
V imply that |I4^(r)| < C 2 \r \ + Cq, we only need to prove that for any t > 0: 


f 1 

lim lim / — W |rj| /" dvR = 0 
£—> 0 n —>-00 / n ^^ 

1=2 


(5.11) 


By the entropy inequality we have: 

[ne] 


kil /r dufs. < ^log + 

i=2 i=2 


Hnit) 


an 


[ne] 


< — y' (^(«) Pi) + G{—oi, Pi)— 20(0, Pi)) + 
na ' 


i=2 


C 

a 


Since Q{a,Pi) + G{—a,Pi) — 2Q{0,Pi) < C'a^, for a constant C independent on i, we 
have 

[ne] 


/ “ Z] I"'*! 

i=2 


C 

< C"ea + —, 
a 


and by choosing a = e (|5.11l) follows. 


□ 
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We are only left to show that 


lim lim 

£—>0 n—>-oo 



[n(l-£)] 


yn 


i=[ne]+l 


^ (-V;V„)G ( - j (y'(r.+i) - (5i)) 


/” du. ds = 0 
(5.12) 


Thanks to the bound (I4.4p . we are now in the same position as in the proof of the 
over-damped dynamics, as considered in [6], and by using similar argument as used 
there (the so called one-block/two blocks) (j5.12ll follows. A slight difference is due to 
the dependence of r on /3j, but since this changes very slowly and smoothly in space it 
is easy to consider microscopic blocks of size k with constant temperature inside. 

At this point the proof of theorem I 2.11 follows by standard arguments. Let Qn 
the probability distribution of tt"’ on C([0,T], A4([0,1]), where A4([0,1]) are the signed 
measures on [0,1]. In appendix B we prove that the sequence Qn is compact. Then, 
by the above results any limit point Q of Qn is concentrated on absolutely continuous 
measures with densities r{y,t) such that for any 0 <t <T, 


f G{y){r{y,t)-r{y,0))dy 
Jo 


-f ^ f ds f G''{y)T{r{y,s),j3{y))dy - G'{1 )t{s) 
Jo Jo 


(5.13) 


= 0 


Furthermore in appendix A we prove that Q is concentrated on densities that satisfy 
the regularity condition to have uniqueness of the solution of the equation. 


6. Appendix A: Proof of the regularity bound \ TW \ 

Proposition 6.1. There exists a finite constant G such that for any limit point distri¬ 
bution Q we have the bound: 


ds J dx {dxT{r{s,x), <C. 


( 6 . 1 ) 


Proof. It is enough to prove that for any function F G C^([0,1]) such that F{0) = 0 the 
following inequality holds: 


ds 


.1 X 1/2 

<C( I F(xfdx] . (6.2) 


dxF'{x)T{r{s,x),fi{x)) — F{1 )t{s) 

In fact by a duality argument, since T(f(s, 1),/3(1)) = r(s), we have: 

^ a, ^^^2 Jo dxF'{x)T{r{s,x),/3{x)) - F{l)f{s) 

/ dx{dxT{r{s,x),/3{x))) = sup - ^ , -• 

Jo FGCUfO.lli L FixYdx 
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Observe that (16.2p corresponds to a choice of test functions G{x) in (I2.19p such that 
G' = F. In order to obtain ()6.2p . compute 

^ n n 

—^n X] + iqi) = ^ F{iln)Alpi 

^ i=l i=l 

n—1 

= ^ Fii/n) {V'in+i) - V'{r,)) + F{1) (f - V\rn)) 

i=l 

1 " 

= - y^VlF{i/n)V'{n) + F{l)f - F{l/n)V'{ri) 

^ i =2 

and after time integration and averaging over trajectories we have 

1 f 


i=l 


/■* /"I /■* 

= / ^^ / + F{1) / f{s)ds (6-3) 

Jo J ^ ^_ 2 J 0 

-F(l/n) l^s I V'in)fyu 0 .. 

It is easy to see that, since F{0) =0 and differentiable, the last term of the right hand 
side is negligible as n —)• oo, by the same argument used in (15.Sp . 

About the first term on the RHS of (|6.3p . by the results of isection 51 it converges, 
through subsequences, to 

_]EQ ii: ds J dxF'{x)T{r{s,x), f3{x))^ . 

About the left hand side of (16.3p . one can see easily that 

1 r ” 

— / - fQ )dvp. ^ 0 . 


2=1 

Using the inequality '^iqf < we can bound the other term of the LHS of 

()6.3p by observing that, for s = 0,t, 




2=1 


\ n ^ 


1/2 


1 2 
ry-t ^ ' 


2=1 


n ^^ n 

2=1 


1/2 


1/2 




1/2 


1/2 


2=1 


2=1 


2=1 


Since U is a continuous function on [0,1] the rhs of the above expression is bounded in 
n and converges to the norm of F as n —>■ oo. Thus (j6.2|) follows. □ 

7. Appendix B: Compactness 

We prove in this section that the sequence of probability distributions Qn on C([0, t], A4) 
induced by 7r„ is tight. Here M. is the space of the signed measures on [0,1] endowed 
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by the weak convergence topology. This tightness is consequence of the following state¬ 
ment. 

Proposition 7.1. For any function G G C^([0,1]) such that G(l) = 0, G'(0) = 0 and 
any e > Q we have 


lim lim sup 

n—>-oo 


sup |< TTnit), G > — < TTnis), G >\ > 8 

0<s<t<T,\s—t\<S 


= 0 (7.1) 

Proof. By doing similar calculations as done in [section 5) (see (15.2p and following ones) 
< TTn{t), G > - < TTn{s), G >= - 

rt 1 n-l 


VnG ( - ) Pi{u) 


/ t •J' ^ 

duY, 

i=i 

‘'V (-) r'(r,+.(u)) - du- 


' V„G () r'(r„(u)) 


n—1 


-h [ du-VnGl-] V'{ri{u)) + -^'^VnG I -] (pift) - Pi{s)) 
Jg 7 V^/ 7^ 


i=2 


^ n _ / . \ 


1 

n ' V j \ n, 

i 

:= /l(s, t) + l2{s, t) + /3(s, t) + /4('S, t) + 15 ( 5 , t) 

We treat the corresponding 5 terms separately. The term I 3 = f* du^VnG (^) V'{ri (u)) 
is the easiest to estimate, since G^(0) = 0, and using Schwarz inequality we have 

G 8 

sup |-^ 3 (s,i)| < sup — / |P'(ri(n))|dn 

0 <s<t<T,|s—1|<5 0<s<t<T,\s—t\<5 Js 

< sup —\t — s\^^‘^( f |P'(ri(M))|^du^ 

0<s<t<T,\s-t\<5 n'J \Js J 

/ 7 AI /2 / rT \ 1/2 

< - ( / \V'{r\{u))\^du 


ni \Jo 


Since, by entropy inequality, 

rT 


E 


'0 


< 


I W(ri(ii))f dn 


1 / 2 - 


< 


Uo 


E(|W(ri(u))p) du 


1/2 


'T/n N 

^ I E [ + 1) 


1 1/2 


du 




so that 


E 


sup 173 ( 5 ,/) I 

0<s<t<T,\s—t\<S 


C'di/2ri/2 

^-172-^ 0- 

777 ,-*-/^ n^oo 


About I 2 , this is equal to 


- iv.G ( 

7 \ n 


du (y'{rn{u)) - r(n)) - duT{u) (7.2) 
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The second term of the above expression is trivially bounded by C5 since |t — s| < <5. 
For the first term on the right hand side of (17.2p . by (1521), we have 


du {V'{rn{u)) - t{u)) 


Pn{s) -Pn{t) 




7 


Pn{u)du + 




n 


{Wn(t) - Wn(s)) 


The last term of the right hand side of the above is estimated by the standard modulus 
of continuity of the Wiener process Wn- For the second term of the right hand side, 
this is bounded by 




rt 


■/ pT \ 1/2' 

E 

sup 7 

/ Pn{u)du 

< 7(5^/2e 

/ Pl{u)du) 


0<s<t<T,|s—1|<(5 

J S 


\Jo J 


< 7^'/' 


L^o 


^(Pn(u))du 


1/2 




Plfudt^pdu 


= ^6^/^ 

1/2 / .T 
0 


HpHu) - Pn^)du + T/3, 


L^o 


3-1 


^ JPndp„fudvp,du + T/d^^ 


{dpnf, 


n\2 

U y 


1/2 


diygduj +T/3, 


-1 


1/2 


1/2 


1/2 




where the last inequality is justified by the inequalities: 




10 


plfyv. < Cn 


" l^^d..du<^ 
~ n 


To deal with the first term we have to prove that 


1 


lim E sup -^\pn{t)\ = 0 

rn-oo \o<t<T?^ / 


(7.3) 


Since 


^ [f{u) - V'irniu))] du 


(7.4) 


The stochastic integral is easy to estimate by Doob’s inequality: 


E sup 

\ 0<t<T 


rt 


27/3T^- / 

n Jo 


CT 
< — 


About the second term, by Schwarz inequality we have that 
E sup 


0<t<T 


f _ V'{rn{u))] du 

Jo 


< 


1 




E (^[r(u) - V'{rn{u))Y^ du 


1/2 
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and by the entropy bound we have 

E([f(u)-y'(rn(u))]') <Cn 

so that this term goes to zero like The first term in (17.4p is trivial to estimate. 

This conclude the estimate of l 2 - 

The estimation of /4 is similar to the proof of (|7.3p . but require a little extra work. 
We need to prove that 


lim E sup 
n->-oo o<t<T 


1 




n—1 





= 0 . 


By the evolution equations we have 


n—l 




i=2 


n—l 


^ ( re ) ~ re2 ^ ( re ) 


i=2 

n-1 




+ / ds ^ - V'{n{s)) 

Jo \nj 

+ l^ds (^V„G(1) W(r„(s)) - V„G V'{r2{s))^ 


(7.5) 


and all these terms can be estimated as in the proof of (17.3p . so that p7.5p follows. 
Also Is can be easily estimated by Doob inequality and using the independence of 

Finally estimating Ji, notice that since G is a smooth function, it can be bounded 
by 


Q 

sup |/l(s,t)| < — sup / 

0<s<t<r,|s—1|<(5 0<s<t<T,\s—t\<5 Js 


< 




7 


r-T >^-1 

/o n 


t n—l 

dre^|W(ri+i(re))| 

i=2 

\ 1/2 


(7.6) 


|l/'(ri+i(re))|^dre j 


and, by entropy inequality 


E 


r-T Y »^-l 


\ 1/2' 

i^*|W(ri+i(re))|2dre I 
'O ™ i=2 / 


< 


i=2 


rT 1 ri-1 




1 " 

>0 ^ i=2 


1/2 




SO that the expression in (|7.6I) is negligible after (5 —>• 0. 


□ 


8. Appendix C: Uniqueness of weak solutions 

Proposition 8.1. The weak solution of (I2.19h is unique in the class of function such 
that 



{dxT{r{x, s), /3{x)))‘^ dx < +oo 


( 8 . 1 ) 
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Proof. Let g{x) > 0 a smooth function with compact support contained in [—1/4,1/4] 
such that J^g{y)dy = 1. Then for A > 0 large enough, define the function 


Gx{y,x) = l- Xg{X{z-x))dz 

J —CXD 

Then for 1/(4A) < a: < 1 — 1/(4A), we have Ga(0,x) = 0 and dyG\{l,x) = 0, and it 
can be used as test function in (I2.19p . So if r{x,t) is a solution in the given class, we 
have 

[ Gx{y,x){r{y,t)-ro{y))dx = -f~^ [ ds f Xg{X{y - x))dyT{r{y, s), I3{y))dy 
Jo Jo Lao 

Letting A —>■ +oo we obtain: 

pX pt 

/ iriy,t)-ro{y))dx = 'y~^ dsdyT{r{x, s),/3{x)), VxG(0,1). 

Jo Jo 

Let ri{x,t),r 2 {x,t) two solutions in the class considered, and define 


R 


px 

{x,t) = / rj{y,t)dy, j 

Jo 


= 1 , 2 . 


By the approximation argument done at the beginning of the proof, we have that 

dtRj{x,t) = -f~^d^T{rj{x,s),/3{x)) 
for every x E (0,1) and t > 0. 

Since r(rj(l,t),/3(l)) = f(t), and since T{r,f3) is a strictly increasing function of r, 
d 

— J {Ri{x,t) - R 2 {x,t)f dx 

= 27 “^ [ {Ri{x,t) - R 2 {x,t))d^{T{ri{x,t),P{x))-T{r 2 {x,t),P{x)))dx 

Jo 

= - 27 "^ / (ri(x,t) - r 2 {x,t)) (r(ri(x,t),/3(x)) - T(r 2 (x, t),/3(x))) dx < 0 . 

Jo 


□ 


9. Appendix D: proof of the entropic hypocoercive bound (4.16) 


We will prove in this appendix that there exists constants A > 0 and C > 0 indepen¬ 
dent of n such that 

j^In{f)<-Xn^In{f) + Gn. (9.1) 

We will use the following commutation relations: 

[dpi,/3^^d*.] = 6ij, [dpi.Al] = = -{d^V dp)i (9.2) 

where d‘^V is the corresponding hessian matrix of V = 

Denote gt = x/W and observe that 

Inigf) = 4 y + \dqgt\% + 2 dqgtQdpgt'j dv^. (9.3) 

Recall that 


in? A 


T,* 

n 




( 9 . 4 ) 
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where 


n—1 


Bn = -n'^ Vnli{i/n)piV'{ri+i) + r?f3{l)pr, 


2=1 


Consequently gt solves the equation: 

dtg = -ri^Al^^'>gt + n^lSngt + i 

9t 2 

We then compute the time derivative of Inidt) by considering the three terms sepa¬ 
rately. The hrst one gives: 


'9t 


A 

dt 


\dpgt\% dv0. =-2n^ j dpgtQdp{A'^^^'^gt) di^p. 


- 2n^7 / dpgtQdp{dp © dpgt) dvfj. 


+ 2n^7 


j dpgtQd^ 


f\dp9t\l 


■^P \ 

V 9t 

+ [ dpgtQdpiBl^^'^gt) dufi.. 


(9.5) 


dvji. 


By the commutation relations (19.21) . and using (|9.4I1 . the first term on the RHS of (|9.5p 
is equal to 

-2n^y dpgtQdggt dv^. -j dp9tQA'"^^^dpgt dv^. 

= -2n^y dpgtQdqgt dvg. - J dpgt®dpgt du^. 

Then the RHS of (19.5p is equal to 

-2n^ / dpgtQdqgt dvp. - 2n^7 J dpgtQdpid* 0 dpgt) du^. 

+ 2n^'y y dpgtQdp dvp. + j gtdpgtQdpB'^^^'> dvp.. 

The last term of the above equation is equal to 

„ „ n—1 

J gtdpgtQdpBn dvp. = -n j gt'^ f3~^VnPi^)V'{ri+l)dp^gt dup. (9.6) 


2=1 


Notice that the term involving n^rpn does not appear in the above expression, because 
the particular definition of ©. For any ai > 0, using Schwarz inequality, (I4.10h and 
()4.11l) . ()9.6I) is bounded by 


J_ /■^ 2 ^'(Vn/ 3 (^))^.,, , 2 , , «in 2 

2q;i 


2=1 


A 


-V'{ri+{f dipp. + 


\^p9t\ps dup. 


Cn ain^ f 2 . 
<— + ^J \dpgt\^ dup. 


for a constant C depending on /3. and the initial entropy, but independent of n. 
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Computing the second term of the RHS of (19.5p we have: 


n 7T. 1 p 

J dpgtQdpid* Q dpQt) dug. = J Yll3]~^\dpdp.g\l dug. + J |9p5|| dug. 

„ n n—l „ 

= / dug. + J \dpg\l dug. 

i=l j=l 


About the third term on the RHS: 

^n—1 Q—1Q—1 

9t 

Summing all together we have obtained 


o „ A.O (\^p9t\%\ _ ^ dp.gt ^p^gt ^p^^p.gt \dpg\Q\dpg\^ 

OpQtQOp I I — -gj2 


9t 


A 

dt 


\dpgt% dpg- ^ j ^pdtOdggt dg^. - r? ^27 - j |(5p3fl| dgf;. 

Cn 


. n—l n 


- 271 % / ^ ^ / 3 j ^ {dp^dp^gt - g^ ^dp^gtdp.gt) dug. + 


j=i i=i 

Now we deal with the derivative of the second term: 


ai 


^ J\dg 9 t\% dug. = j dqgtQdqiA^^^^gt) dug. - 2n^7 J dqgtQdq{d* Q dpgt) dug. 

+ 271^7 / dqgtQdq 


) dug. + / dqgtQdq{Bngt) dug. 


. n—l n 


/ nit - X It 

dqgtQdq{A"^*^ gt) dug. - 277 % J {dpidg^g - gi^dp.gdq.g)'^ dug. 


j=l i=l 


+ J dqgtQdq{Bngt) dug.. 

(9.8) 


The first and the last term give: 

- 277 ^ y dqgtQdq{A^^^'>gt) dug. + J dqgtQdq{Bngt) dug. 

= 2n^ j dqgtQ{dgVdp)gt dug. + j gtdqgtQdqBn dug. 
The last term on the RHS of the above expression is equal to 
j gtdqgtQdqBn dug. 


n—l 


^dtidq^gt) 


i=2 


* '' \„ •ri a ( 1 


77 


V,7/l - V"{ri+i)pi - VnP 


77 


V (ri)pi_i 


dua 


+ n J13.^ V"{r2)pi. 
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Since V" and V„/3 are bounded and j3{-) is positive bounded away from 0, this last 
quantity is bounded for any 02 > 0 by 

/ « n—1 « 

\dqgt\%d-iyi3. + Ca2^ / '^pfgfdvp. <n^a2 / lOggil^^du/s. + C'a^^n. 
i=i 


Since V" is bounded, for any 0:3 > 0 we have 

2n^ j dqgtQ{dlVdp)gt dvj^. < J |(9gC/tl| diyp. 


+ 


I ^ loo'*' 

as 


j \dp9t\Q dvp. 


Putting all the terms together, the time derivative of the second term is bounded by 

Cn^ 


A. 

dt 


J \dqgt\% dup. < {a2 + a'i)r? J Idqgtf^ dv0. + J |Sp5t|| dup. 


. n—1 n 


2n^7 ^ {dp^dq.g - g^ ^dp^gdq.gf dup. + C'a^ 

J ^ _ 1 _T 


(9.9) 


j=ii=i 


About the derivative of the third term, using the third of the commutation relations 
([ 0 |) . gives 


dt' 


dvp. 


' I dqgtQdpgtduf,. =-2n^ /I dq{A"^^'> gt)Qdpgt + dqgtQdp{A^^^'> gt) 

j \dq(yBYigt)Qdpgt “b dqgtQdpi^Bng^')^ dup. 

-2n^7 f [dqgtQdp{dp © dpgt) + dqid* 0 dpgt)Qdpgt] dup. 


+2n^7 


( \dpgt\%i \ ^ Q ( I dpQt 11 

" V 9t 


\ 9t 


dqQt&dp 

= 2r? jidqVdp)gtQdpgtdi^p. - 2n‘ 
1 


Qdpgt 


dvp. 


J l^ggil^dl'p. 


“1“ 2 y 9t [pqBfiQdpgt -h dqgtQdpB’n^ dvp. 

« n—1 n 

-dnS [{dpAjd) {dpAj9)] 

j=l i=l 

P n—1 n 

+2n^7 I ^ [idpAj9)idpi9t)idq,gt) + idqAi9)idpj9t)idp,9t)] dvp. 

j=i i=i 

P n—1 n 

—4n^7 I ^^J^9 Hdp,9tfidq.gt)idp.gt)dup. 

j=l 2=1 

(9.10) 


The last three terms of the RHS of the (I9.10p can be written as 

/ n—1 n 

[{AA9t - gi^dp^gtdq.gt) {dpAj9t - gi^dp^gtdp.gt)] dup. 

j=i i=i 
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SO they combine with the corresponding terms coming from the time derivative of the 
first two terms of In giving an exact square. 

The second term of (I9.10p . by the same arguments used before, can be bounded by 


n 


'a^ J \dqgt\'^^dui 3 . + j \dpgt\\di'p. + Cn{a^ ^ + 0 : 5 ^) 


About the first term of (|9.10l) . since V" is bounded, it is bounded by f \dpgt\'^^dvp.. 
Putting all these bounds together we obtain that 

^In{ft)<-n^Kp j Idpgtf^dup. - Kq j \dqgt\'^^diyp. - 2n^ j dpgtQdqgtdu^. + Cn 


. n—1 n 


-2N^-f / ^ ^/ 3 . ^ [{dpidq.gt - g^ ^dp.gtdq.gt) + {dp.dp^gt - gt ^dp.gtdp^gt)] dvfi. 

j=i i=i 


with 


«i c 

Kp = 2-/ - — -as - Vc 

as 


1 / 

00 


-P 2 

Kq = 2 — Oi2 — as — 04 

By choosing 02 + as + 014 < 1 we have obtained that for some constants Ci, C 2 > 0 
independent of n 

^Inift) <-n^In{ft) + Cin + C2 f Idpgtll^duf}.. 


By recalling that 

J j 

after time integration we have for some constant C 3 : 


, c' 

dvp. < — 
^ n 


that implies 


/n(/i)<e-"'‘4(/o) + —(1-e-"'') 

n 


Uft) < — 
n 


(9.11) 


for any reasonable initial conditions such that In{fo) is finite and not growing too fast 
with n. 

Remark 9.1. An important example for understanding the meaning of a density with 
small In functional, eonsider the inhomogeneous Gibbs density: 

( n n—1 ^ \ 

^ Pinri + ^ -yn{PiTi)pi j /M (9.12) 

j=i i=i ^ ) 

where M is a normalization constant. In the case of constant temperature these densities 

play an important role in the relative entropy method (cf [141 15] ). as to a non-constant 
profile of tension corresponds a profile of small damped velocities averages. Computing 
In on f we have 

n—1 


InU) = E 


2=1 


fdiTi /^i+lTj+l T Vn(/3iTi) 
n 


= 0 . 
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